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ABSTRACT: We present our analytic results for the NLO corrections to the partial decay 
width t —> H + b followed by b — > BX for nonzero b-quark mass (m& 7^ 0) in the Fixed- 
Flavor-Number scheme (FFNs). To make our predictions for the energy distribution of 
the outgoing bottom-flavored hadron (B-hadron) as a function of the normalized B-energy 
fraction xb, we apply the General-Mass Variable-Flavor-Number scheme (GM-VFNs) in 
the general two-Higgs-doublet model. In order to describe both the b-quark and the gluon 
hadronizations in top decay we use fragmentation functions extracted from data from e + e~ 
machines. We find that the most reliable prediction for the B-hadron energy spectrum is 
made in the GM-VFN scheme. 
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1 Introduction 

The top quark has been the latest standard model particle discovered by the CDF and DO 
experiments at Fermilab Tevatron [1]. Its properties are essential for our understanding of 
the standard model (SM) theory. At the Large Hadron Collider (LHC), one expects a cross 
section a(pp — >■ tiX) ~ 1 (nb) at design energy \fS = 14 TeV [2]. With the LHC design 
luminosity of 10 34 (cm) _2 (sec) _1 , it is expected to produce a tt pair per second. Thus, the 
LHC is a superlative top factory, which allows to carry out precision tests of the SM and, 
specifically a precise measurement of the top quark properties such as its mass. Due to 
the element \Vtb\ ~ 1 of the Cabibbo-Kobayashi-Maskawa (CKM) [3] quark mixing matrix, 
the top quark is decaying dominantly through the mode t —> bW + within the SM. As it 
is well known, bottom quarks hadronize before they decay, therefore each b-jet contains 
a B-hadron which most of the times is a B-meson. Events with B-mesons are identified 
by a displaced decay vertex associated which charged-lepton tracks. This is precisely the 
signature used to identify b-jets. In Ref. [4], we studied the B-meson energy distribution 
produced in top decay considering both the bottom quark and the gluon fragmentations. 
We also studied the angular distribution of the W-boson decay products in the decay chain 
t bW + — > Bl + ui + X. The effects of b-quark and hadron masses are also considered. 

In many extensions of the SM such as the minimal supersymmetric standard model 
(MSSM), the Higgs sector of the SM is enlarged, typically by adding an extra doublet of 
complex Higgs fields. After spontaneous symmetry breaking, the two scalar Higgs doublets 
Hi and Hi yield three physical neutral Higgs bosons (h, H, A) and a pair of charged- 
Higgs bosons (H^) [5]. The top quark is decaying dominantly through t — > bH + in a 
model with two-Higgs-doublet (2HDM) [6], providing that the top quark mass mt, bottom 
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quark mass m& and the charged-Higgs boson mass mjj satisfy m< > m& + rntr. In this 
case one expects measurable effects in the top quark decay width and decay distributions 
due to the i/ ± -propagator contributions, which are potentially large in the decay chain 
t -> bH + -> b(r + v T ). 

At LHC, the dominant source of top quarks is pp — > ti process, therefore the charged Higgs 
boson has been searched for in the subsequent decay products of the top pairs 
and ti — > H^H^bb when decays into r lepton and neutrino. 

In our previous work [7], we studied the energy spectrum of the inclusive bottom-flavored 
mesons in the presence of charged Higgs boson by working in the massless scheme or zero- 
mass variable-flavor-number (ZM-VFN) scheme where the mass of bottom quark is set to 
zero from the beginning. As it was shown, in the limit of vanishing b-quark mass our results 
in two variants of the 2HDM are the same. In the present work, we impose the effect of 
b-quark mass on the B-spectrum employing the general-mass variable-flavor-number (GM- 
VFN) scheme. As it is shown, the results will be different in two variants of the 2HDM and 
it is found the NLO corrections with mj 7^ to be significant. 

To obtain the total B-hadron energy distribution of the top decay, two contributions due 
to the decay modes t — > bW + (in the SM) and t — > bH + (in the 2HDM) should be summed 
up. However, the contribution of SM is always larger than the one coming from 2HDM [7], 
but there is a clear separation between the decay channels t — > bW + and t — > bH + in both 
the tiX pair production and the t/iX single top production at the LHC [8]. New results of 
a study on the charged-Higgs bosons in pp collision at a center of mass energy of y/s = 7 
TeV are reported by the ATLAS Collaboration [9]. 

Finally, since bottom quarks fragment into the B-meson, in order to describe the fa- 
quark non-perturbative fragmentation some phenomenological hadronization models can be 
used. Our treatment at NLO in the GM-VFN scheme is manifestly based on the factoriza- 
tion theorem of QCD which guarantees that the fragmentation functions are universal and 
subject to DGLAP evolutions [10]. Relying on the universality of the hadronization mech- 
anism, we can tune such models to data on B production in e~e + annihilation data from 
CERN LEP1 and SLAC SLC and use them to predict the B-hadron spectrum in top decay. 
The hadronization of the b-quark was considered in the NLO QCD analysis of top-quark 
decay in Refs. [11, 12] and was identified to be the most important factor of uncertainty in 
the determination of the top-quark mass. 

This paper is organized as follows. In Sec. 2, we give the parton-level expressions for the 
NLO QCD corrections to the tree-level rate of t — > bH + in the fixed flavor number (FFN) 
scheme. In Sec. 3, the scheme of GM-VFN is explained by introducing the perturbative 
fragmentation function 6—7-6. In Sec. 4, we present our hadron-level results working in the 
GM-VFN scheme. In Sec. 5, we summarize our conclusions. 

2 Parton level results 

2.1 Born level rate of t — > bH + 

We consider the decay channel t — > bH + in the general 2HDM, where H\ and H2 are the 
doublets whose vacuum expectation values give masses to the down and up type quarks, 
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respectively, and a linear combination of the charged components of Hi and H2 gives the 
physical charged Higgs H + . In a general model with two Higgs doublets to avoid tree level 
flavor-changing neutral currents (FCNC), the generic Higgs coupling to all quarks has to be 
restricted. There are two possibilities for the two Higgs doublets to couple to the fermions. 
In the first possibility (model I), one of the Higgs doublets (Hi) couples to all bosons and 
the other one (#2) couples to all the quarks. In this model, the Yukawa couplings between 
the charged Higgs boson, the top and the bottom quarks are given by [13] 



Li = 2 J^ mw V tb p\H + t [mt(l - 7s) - m 6 (l + 75)] &} 



+ H.c. (2.1) 



For the vacuum expectation values of Hi(v\) and -#2(^2)1 we have v 2 + v 2 = (\^2G F )~ 1 
where G F is the Fermi's constant and the ratio of the two values is a free parameter and 
one can define the angle (3 to parameterize it, i.e. tan/3 = V2/V1. The weak coupling factor 
gw is related to the Fermi coupling constant by g 2 ^ = Ayplm^Gp. In the equation above, 
H+ = cos fiHy — sin fiHf is the physical charged Higgs boson. 

In the second possibility (model II), the doublet H2 couples to the right-handed up-type 
quarks (un,CR,tn) and the Hi couples to the right-handed down-type quarks. In this 
model, the interaction Lagrangian would be 

L11 = ^— V tb ^H+t[mt cot (3(1 - 75 ) + m b tan 0(1 + 75)] 6 J +H.c. (2.2) 

The Born amplitude for the process t — > bH + can be parameterized as Mq = u&(a + 675)1^, 
therefore, the tree-level total decay width is given by 

f„ = + 4>)[1 + 4 - + «(•* - W, 4. (2-3) 

1D7T [ mf mf mt J mf mf 

\2 



where X(x, y,z) = (x — y — z) — Ayz is the Kallen function and for model /, one has 

a 2 + b 2 = V2\V tb \ 2 G F (m 2 + m 2 b ) cot 2 /3, 



a 2 



b 2 = -2y/2\V th \ 2 G F (m h m t ) cot 2 /3, (2.4) 



and for model II 



a 2 + 6 2 = V2\V tb \ 2 G F (m 2 cot 2 ft + m 2 b tan 2 ft), 

a 2 -b 2 = 2V2\V tb \ 2 G F (m b m t ). (2.5) 



In the limit of vanishing b-quark mass, the tree level decay width is discussed in [7]. Since 
m b <C mt, the b-quark mass can always be safely neglected in model I but in model II, 
the left-chiral coupling term proportional to m b tan /3 can become comparable to the right- 
chiral coupling term mt cot (3 when tan (3 becomes large. Therefore, one can not naively set 
m b = in all expressions in model II. 

In the following, we explain the calculation of the NLO QCD corrections to the Born 
level decay rate of t — > bH + and we present the parton-level expressions for dT(t — > BH + + 
X)/dxB at NLO in the FFN scheme where m b ^ is considered. 
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2.2 Virtual one-loop corrections and counterterms 

The virtual corrections to the tbH + -veitex consists of both infrared (IR) and ultravi- 
olet (UV) singularities where the ir- and uv-divergences arise from the collinear- and 
the soft-gluon singularities, respectively. In our calculation, we adopt the on-shell mass- 
renormalization scheme and all singularities are regularized by dimensional regularization 
in D = 4 — 2e space-time dimensions. To simplify the formulas we introduce the follow- 
ing kinematic variables, in the notation of Refs. [14, 15], along with some other required 
variables 

po = -(1 + R-y), 



Po 



Ps = PoV 1 ~ P 2 , 

P± = PO ±P3, 



v 1 l P+ 
2 p_ 

Y w = - In ■ 



2 l-p+ 
H = (a 2 + b 2 ) Po + (a 2 - b 2 )VR, 



T = Po(l — Xb)\/xl — f3 2 , 



Hx b ) = po [y/x 2 b - p 2 - In J—] . ( 2 - 6 ) 



°°b ~ \] x l ~ P 2 

where the scaled masses R = m 2 /m 2 and y = m 2 H+ /m 2 are defined. Choosing these 
notations, the tree-level total width (Eq. (2.3)) is simplified to Tq = mtHp^/(4Tr). It is 
also convenient to introduce the normalized energy fractions Xj = Ei/E™ ax {i = b,g) where 

jpmax 

t, h — m t po- 

Considering the above notations, the contribution of virtual corrections into the differential 
decay width reads 

^vir 



drr p 3 



dXh 8-7TTO/ 



\M vir \ 2 5(l- x b ), (2.7) 



where, |M vir | 2 = l/2J2 Sp in( M o M loo P + M i oop M o)- The renormalized amplitude of the 
virtual corrections is written as Mi oop = Ub(A c t + A;)ut, where A C f stands for the counter 
term and A; arises from the one-loop vertex correction. Following Refs. [14-16], the counter 
term of the vertex includes the mass and the wave-function renormalizations of both the 
top and bottom quarks as 

A ct = (a + b)(^ + ^-^Y^ + (a-b)(^ + ^-^Y-^, (2.8) 
v ' \ 2 2 m t J 2 v ' \ 2 2 m b J 2 ' v ' 
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where, the mass and the wave function renormalization constants read 
5m q _ a s (//fl) i 



5Z n 



4?r 

47T 



C F 



— -3 7£ ; + 31n^^ + 4 



-^ + -^-3 7i?+ 31n^# 



+ 4 . 



(2.9) 



Here, m q is the mass of the relevant quark and C F = (N 2 -l)/(2N c ) = 4/3 for N c = 3 quark 
colors. In the equation above, em and ejjv represent infrared and ultraviolet singularities 
and "fE = 0.577216 • • • stands for the Euler constant. 
The real part of the one-loop vertex correction reads 



A, 



7T 



C F [(a 2 - b 2 )VRG + + (a 2 + b 2 )G. 



(2.10) 



with 

G+ = ^m 2 p C {ml,m 2 ,?n 2 H+ ,ml,0,m 2 ) + 5 (m 2 , 0, m 2 ) + 2B (m 2 H+ ,m 2 ,m 2 ) 

+B (m 2 ,0,m 2 )-2, 
G- = Am 2 plC (ml,m 2 ,m 2 H +,ml,0,m 2 ) + (2p - R)B (ml,0,ml) 

+(1 + R)B (m 2 H+ ,m 2 b ,m 2 ) + (2p - l)B (m 2 , 0, m 2 ) - 2 Po , (2.11) 

where, Bq and Cq functions are the Passarino-Veltman 2-point and 3-point integrals which 
can be found in Rcf. [17]. 

All uv-singularities are canceled after summing all virtual corrections up but the ir-divergences 
are remaining which are now labeled by e. Putting everything together, the virtual differ- 
ential decay rate normalized to the Born total width, reads 



1 (IT™ 



f dx b 2tt 



-C F 5(l-x b ){ -2 + hiR 



2(1 -Po) _ 3abp _ 2po v 

y H p 3 



P3 



Y p K-lny 



ypoH 



[2^R{a 2 - b 2 ) + {a 2 + 6 2 )(1 + R)} \ 



P3 



111 



47T/4 



, 1 \ 2 P0 

IE + 



Li 2 (p-) ~ Li 2 {p+) + Li 2 (l ) 

P+ 



(2.12) 



where, Li2(x) = — J^(dt/t) ln(l — i) is the Spence function. In Eq. (2.12), one has 

Gp 

ab = -j=\V t b\ 2 {m 2 t - m 2 ) cot 2 /3, 



(2.13) 



in model I, and 



Gf 

ab = — — | V t b\ 2 (m 2 cot 2 /3 — m 2 tan 2 /3), 
v2 



(2.14) 



in model II. The terms a 2 + b 2 and a 2 — b 2 are given in Eqs. (2.4) and (2.5) in both models. 
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2.3 Real one-loop corrections (Bremsstrahlung) 

The 0(as) real gluon emission amplitude reads 

M^^gs^u^s,) ^-^ - 2P 1 + ^ 9 }(a + b l5 )u(p t ,s t )e;(p g ,r), 

A "Pi ' Pg ^Pb ' Pg 



(2.15) 



where the polarization vector of the gluon is denoted by e(p g ,r). As before, to regulate the 
IR-divergences we work in D = 4 — 2e dimensions and for simplicity, we choose the top 
quark rest-frame. To get the correct finite terms in the normalized differential decay rate, 
the Born width To will have to be evaluated in the dimensional regularization at 0(e), i.e. 
To ~~ * To{l — e[21n(2^3) + 7^ — ln(47r/i 2 /m 2 ) — 2]}. When one integrates over the phase 
space for the real-gluon radiation, terms of the form (1 — x b )~ 1 ~ 2e arise which are due to the 
radiation of a soft gluon in top decay. Therefore, for a massive b quark, where x^ m i n = (3, 
we use the following expansion 

±5(1 -x b ) + - 1 . +0(e), (2.16) 



(l-s 6 )i+* 2e " °> ' (1-3*)+ 
where the plus distribution is defined as 

1 f(Xb) dXh =t f{ - Xh }- m dx h + fil)Hl-P). (2.17) 



2.4 Parton-level results for dT/dx a in FFN scheme 

Now, we present our analytic results for partial decay rate normalized to the Born width 
in the FFN scheme, by summing the tree level, the virtual and the real contributions. Our 
result reads 



1 dT b C F a s (ix R ) j xn , 
- <5(1 - x b ) H { (5(1 - x h ) 



T d,X b 7T 



[1-R)Y W 2 In 2poP+ l 



P3 -Jv 



-2— |L2 2 (p_) -Li 2 (p + ) + Li 2 (l ) +2-^ p 3 H hpohi h 

P3 P+ P3 V 2 P+ 

-4-[2(a 2 - 6?)VS+ (1 + R)(a 2 + 6 2 )]) + ( - + V ~ P ° - ^Y w - ^ab) InR 
2yH L J V P3 2H ' 

_ 2(T + x b *(x>)) [VM 2 _ ft2 + 2 + &2 j _ pQ a*+^ [T+ ^ 1 

•Hp 3 (l - ^f>) + #P3 J 

(2.18) 

Integrating dT b /dx b of Eq. (2.18) over x&(/3 < x b < 1), we obtain the total decay rate 
presented in Refs. [14, 15]. 

Since, the B-meson can be also produced from the fragmentation of the emitted real 
gluon, we also need the differential decay rate dT g /dx g in the FFN scheme, where x g = 
Eg/(mtPo) is the scaled energy fraction of the real gluon. To calculate the dT g /dx g , we 
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integrate over the momentum of b-quark by fixing the momentum of gluon in the phase 
space. Our result is listed here 



1 dV g _ C F a s (n R )pl(l 
f o dx g TtP3XgH 




1-Xg 



1 



) + 2/3(a 2 - b 2 ) 



x tanh 




1 - 



/3 2 (1 - 2p x g ) _ I" 2 2 



2 + 



p x 2 g [2-p {l + 3x g )] 
(1 - 2p x 9 ) 2 



) 



+2/3(a 2 -6 2 ) Jl 




(1-X,)2 J" 



(2.19) 



3 GM-VFN scheme 

Our main purpose is to calculate the scaled-energy distribution of the B-hadron produced in 
the inclusive process t — > BH + + X in the 2HDM, where X stands for the unobserved final 
state. Thus we calculate the partial decay width of the corresponding process differential 
in xb { dT/dxB), at NLO in the GM-VFN scheme, where xb = EB/(mtPo) is the scaled 
energy fraction of the B-hadron. In the top quark rest frame, the B-hadron has energy 
Eb = Pt " PB/ m t, where tub < Eb < [m 2 + rn B — m 2 H+ \/ (2rrit). In the case of gluon 
fragmentation, g — > B, it has energy mg < Eb < [m 2 + m 2 B — {nib + m H+) 2 ]/ {^ m t)- 
According to the factorization theorem of QCD [18], the B-hadron energy distribution 
can be obtained by the convolution of the parton-level spectrum with the nonperturbative 
fragmentation function (z, ^p), 



where, [ip is the factorization scale and hr is the renormalization scale which is related to 
the renormalization of the QCD coupling constant. A choice often made is to set u,r = fj,p 
and we shall use this convention for most of the results. Here, dT^ M /dx a is the differential 
decay width of the parton-level process t — > a + X at NLO in the GM-VFN scheme, where 
X comprising the H + boson and any other parton. We now discuss the evaluation of the 
quantities dF^ M /j,p)/dx a in the GM-VFN scheme in detail. 

In Ref. [7], using the ZM-VFN scheme we evaluated the quantities 1/To X dT a /dx a {a = b, g) 
for the process t — > a + X, where nif, = is put right from the beginning. In this scheme, nib 
only sets the initial scale [i 1 ™ = 0{nib) of the DGLAP evolution, however all information 
on the nib dependence of dT a /dx a is wasted. 

In sec. 2 of the present paper, we applied the FFN scheme which contains of the full nib 
dependence. In this scheme, the large logarithmic singularities of the type (a s /-7r) In R, 
where R = m 2 /m 2 , spoil the convergence of the perturbative expansion when nib / m t 
(see Eq. (2.18)). The GM-VFN scheme is devised to resum the large logarithms in nib an d 
to retain the whole nonlogarithmic nib dependence at the same time. This is achieved by 
introducing appropriate subtraction terms in the NLO FFN expressions for dT a /dx a , so 
that the NLO ZM-VFN results are exactly recovered in the limit mb/nit — > 0. With this 



dT 

dxB 



a=b,g 




{fl R ,U,p)D a ( ,flp) 
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explanation, the subtraction terms are constructed as 



1 dT 



Sub 



1 dt FFN 1 df ZM 
lim a 



Tq dx a m b -s-0 Tq dx a 



T dx a 



(3.1) 



and the GM-VFN results are obtained by subtracting the subtraction terms from the FFN 
ones [19, 20], as 



i dr? M i df FFN i dT? ub 



T dx a T dx a T dx a 
Taking the limit m b — > in Eqs. (2.18) and (2.19), we obtain the subtraction terms as 



(3.2) 



1 dT^ ub 
To dx b 



1 dT s g ub 
To dx g 



2tt 



-C f 



l-x b 



In 



/4 

2 

mt 



21n(l -x 6 ) - 1 



R, 



2tt 



C f 



l + (l-x B 



In — ^ — 2 In x Q — 1 

m 2 



(3.3) 



(3.4) 



As it is guaranteed by Collin's factorization theorem [18], the subtraction terms are universal 
and as we already presented in Ref. [4], Eq. (3.3) coincides with the perturbative FF of 
the transition b — )• b [21]. 



4 Numerical results 

Now we present our phenomenological predictions by performing a numerical analysis. As 
it is referred in Ref. [22], a charged Higgs having a mass in the range &0GeV < mjj± < 
IGOGeV is a logical possibility and its effects should be searched for in the decay modes 
t — > bH + — > Bt + v t + X. A beginning along these lines has already been made at the 
Tevatron [23, 24], but a definitive search of the charged- Higsses over a good part of the 
(rriff+ — tan /3) plane is a program that still has to be carried out and this belongs to the 
LHC experiments [25]. 

Following Ref. [26], we adopt the present lower limit mjj+ > 79.3 GeV obtained from LEP. 
From Ref. [26], we use the input parameter values Gf = 1.16637 x 10~ 5 GeV -2 , mt = 
172.0 GeV, m b = 4.90 GeV, m B = 5.279 GeV, and \V tb \ = 0.999152. We evaluate a^Qm) 
at NLO in the MS scheme using 



( n /) i \ 



with bo and b\ given by 



1 



6 log(/i 2 /A 2 ) 



1 



6 1 log[log(/x 2 /A 2 )] ^ 
6 2 log(^/A 2 ) V 



33 - 2n f 
12^ : 



153 - I9n f 
24^ : 



(4.1) 



(4.2) 
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Figure 1. dT(t —> BH + + X)/dxB as a function of xb in the GM-VFN scheme for model I, taking 
m H+ = 120 GeV and tan /3 = 1, 5 and 10. 

(5) 

where A is the typical QCD scale and we adopt Al^i. = 231.0 MeV adjusted such that 

a s = 0.1184 for mz = 91.1876 GeV [26]. In Eq. (4.1), rif is the number of active quark 
flavors. 

To describe the transitions b,g — > B, we employ the realistic nonperturbative S-hadron 
FFs determined at NLO in the zero-mass scheme through a global fit to e + e _ -annihilation 
data presented by ALEPH [27] and OPAL [28] at CERN LEP1 and by SLD [29] at SLAC 
SLC. Specifically, for the b — )■ B transition the power model D^z^fi 1 ™) = Nz a (l — z)" was 
used at the initial scale f/p 1 = 4.5 GeV, while the light-quark and gluon FFs were generated 
via the DGLAP evolution. The result of fit read N = 4684.1, a = 16.87, and /3 = 2.628 
[30]. Note, if the same experimental data are fitted in the ZM-VFN and GM-VFN schemes, 
the resulting FFs will be somewhat different. As it is shown in Ref. [31], the hadronization 
of the bottom quark is identified to be the largest source of uncertainty in the measurement 
of the top quark mass. 

To present our results for the scaled-energy (xb ) distribution of B-hadrons, we consider 
the quantity dT(t -> BH + + X)/dx B taking the H + boson to be stable. In Fig. 1, we show 
our prediction for the size of the NLO corrections in the GM-VFN scheme for the model 
I. Here, the mass of Higgs boson is fixed to mjj+ = 120 GeV and the different values of 
tan/3 are considered, i.e. tan/3 = 1,5 and 10. However, as in Ref. [32] it is pointed out, 
the small values of tan /3 are excluded by the indirect limits in the (m^±,tan/3) plane. For 
example, taking the CP-conserving scenario m^-max and a top quark mass of 174.3 GeV , 
values of tan /3 between 0.7 and 2.0 are excluded, but this range depends considerably on 
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X B 



Figure 2. dT(t —> BH + + X)/<1xb as a function of xb for model I. Different values of the Higgs 
boson mass is considered, i.e. m#+ = 80, 100 and 120 GeV. Other free parameter is fixed to 
tan f3 = 5. 

the assumed top quark mass and may also depend on Msusy (the soft SUSY breaking 
scale parameter). 

In Fig. 1, Both the b-quark and gluon fragmentations are included. In Ref. [7], we showed 
the g — >• B contribution is negative and appreciable only in the \ow-xb region. For higher 
values oi xb the NLO result is practically exhausted by the b — >■ B contribution. 
From Fig. 1, it can be seen that when tan/3 is increased the decay rate is decreased and 
the peak position is shifted towards higher values of xb- In Ref. [14], it is shown when the 
values of tan/3 exceed tan/3 = 2, the decay rate becomes quite small. Here, the mass of 
B-hadron creates a threshold at xb = 2ms/(mt(l + R — y)) = 0.12. 

Adopting the limit 80GeV < m H ± < 160GeV from Ref. [22], in Fig. 2 we study the 
energy spectrum of the B-hadron for model I in different values of the Higgs boson mass, 
i.e. Tnjj+ = 80,100 and 120 GeV, by fixing tan/3 = 5. As mentioned, the mass of B- 
hadron is responsible for the thresholds at xb = 0.08 (for m H + = 80 GeV), xb = 0.09 (for 
m H + = 100 GeV) and x B = 0.12 (for m H + = 120 GeV). 

In Fig. 3 and Fig. 4, the results presented in Figs. 1 and 2 are regained for model II. 
The thresholds are as before. 

In Fig. 5, we compared the energy spectrum of B-hadron in the GM-VFN (rrif, ^ 0) 
and ZM-VFN [mi, = 0) schemes, using tan/3 = 1 and rriH+ = 120 GeV. As it is shown, 
the results of both models are the same in the GM-VFN scheme but the result of ZM-VFN 
scheme shows an enhancement in the size of decay rate about 1.3% at xg = 0.8. In Figs. 6 
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X B 



Figure 3. xb spectrum in top decay considering the decay mode t — > BH + +X, taking m H + = 120 
GeV and tan/3 = 1,5 and 10 in model II. Threshold at xb is shown. 




X B 



Figure 4. xb spectrum for model II considering the different values of the Higgs boson mass, i.e. 
m H + = 80, 100 and 120 GeV, by fixing tan/3 = 5. 




Figure 6. As in Fig. 5, but using tan/3 = 5. 



and 7, this comparison is done using tan/3 = 5,10 and m H + = 120 GeV. As it is seen, 
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the prediction for the energy spectrum extremely depends on the model at the large values 
of tan/3 when mj, ^ 0. Therefore, our most reliable prediction for the energy spectrum is 
made at NLO in the GM-VFN scheme. 

In Fig. 8, the energy spectrum of B-hadron in decay modes t — > BW + + X and t — > 
BH + +X are compared. As before, the mass of B-hadron create the thresholds at xb = 0.12 
(for m H + = 120 GeV) and x B = 0.08 (for m w + = 80.399 GeV). It is obvious that the 
contribution of the top decay mode in the SM is always larger than the one coming from 
the 2HDM. However, to obtain the total energy spectrum of B-hadron in the top quark 
decay all decay modes including t — > B + W + /H + should be summed up. 



5 Conclusions 

Clearly, the decay modes t —> W + + B have been and will be the prime source of information 
on the top quark mass. We have studied these dominant decay modes along the lines of 
Ref. [4]. In the theories beyond the standard model including the two-Higgs-doublet, the 
top quarks also decay into a charged Higgs and a bottom quark thus it may be useful to 
also use t — > H + + B events for a cross check. The search for the light charged Higgs boson 
(rriff+ < mt) produced from the decay mode t — > bH + (— > t + v) has been performed using 
about lfb~ l of data collected in proton-proton collisions at t/s = 7 TeV [5]. 
In our previous work [7], we applied the ZM-VFN scheme to study the dominant decay 



-13- 



> 4 

a 



X 
T3 



X 
+ 



T3 



NLO-Model I (t -> BH + +X) 
NLO-Model II (t -> BH + +X) 
NLO (t -» BW + +X) 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

X B 

Figure 8. Xb spectrum in top decay considering the decay modes t —> BW + + X (solid line) 
and t — > BH + + X (dashed and dotted lines), taking m w + = 80.399 GeV, m H + = 120 GeV and 
tan/3 = 1. 



channel t — > BH + + X in the 2HDM, where the mass of bottom-quark was set to zero 
and thus all information on the mj dependence of the B-hadron spectrum was wasted. In 
the present work, we studied the energy spectrum of B-hadron in top decay considering 
the quantity dT/dxB in GM-VFN scheme. Our main purpose was to investigate both the 
effect of b-quark mass and the gluon fragmentation to the B-hadron energy distribution. In 
order to study these effects we have calculated, for the first time, an analytic expression for 
the NLO radiative corrections to the differential top decay width dT/dx a (a = b,g) in two 
variants of the 2HDM. To ensure our results, we have checked that by integrating dT/dxi, 
over Xfe we recover the known results presented in Refs. [14, 15]. 

In Ref. [7], we showed in the limit m& — > our results in both models are the same but 
in the present work we have checked the results for the energy spectrum of B-hadron are 
completely different for two models in the general 2HDM when nib ^ and tan (3 is large. 
In conclusion, the most reliable predictions for dT/dxB is made at NLO in the GM-VFN 
scheme. 

Comparing future measurements of dT/dxs at the LHC with the presented predictions will 
be important for our understanding of the Higgs coupling in 2HDM and new physics beyond 
the standard model. Our results in Refs. [4, 7] and the present work both will be able to 
test the universality of the B-hadron fragmentation functions and provide a clean method 
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to gauge the normalization of Monte Carlo event generators. 
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